Generalized Ward identity and gauge invariance of the color-superconducting gap 
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We derive a generalized Ward identity for color-superconducting quark matter via the functional 
integral approach. The identity implies the gauge independence of the color-superconducting gap 
parameter on the quasi-particle mass shell to subleading order in covariant gauge. 
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Due to asymptotic freedom of quantum chromodynamics (QCD) the strong coupling constant g becomes small at 
large momentum Q. In cold, dense quark matter, the typical momentum scale is of the order of the quark chemical 
potential (i. Consequently, at asymptotically large densities, where \i 3> Aqcd, quark matter is a weakly coupled 
system. In this case, the dominant interaction between two quarks is single-gluon exchange, which is attractive in 
the color-antitriplet channel. Due to Cooper's theorem ; 2j, at sufficiently low temperatures the quark Fermi surface 
becomes unstable with respect to the formation of Cooper pairs. Since this is analogous to what happens in ordinary 
superconductors |2|, this phenomenon was termed color superconductivity [3|- Recently, color superconductivity has 
been extensively studied within QCD at weak coup lin g jj l la, IE S j m phenomenological Nambu-Jona-Lasinio-type 
models 0,113, and in effective theories 111! Il2l Il3l Il4l llfil. 

In a superconductor, exciting particle- hole pairs costs at least an energy amount 20Oj where 0o is the value of the 
superconducting gap parameter at the Fermi surface for zero temperature. The gap parameter can be computed from 
a gap equation derived in mean-field approximation which involves single-gluon exchange and bare quark-quark-gluon 
vertices. Schematically, this gap equation can be written in the form 



4>a = g 

In weak coupling, g -C 1, the solution is 
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6 = 2 6/iexp( -- ) [1 + 0(3)] 
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The first term in Eq. contains two powers of the logarithm ln(/x/(/) ): one is the same as in BCS theory (the 
so-called "BCS logarithm") and the other arises from the collinear exchange of almost static magnetic gluons, which is 
a long-range interaction |3, |6( . The weak-coupling solution J5| implies that this term contributes to the gap equation 
at order 0(1). We call this term the leading-order term. The value of the coefficient £ determines the constant c in Eq. 
|2"|l. The second term in Eq. contains subleading contributions of order 0(g) to the gap equation, characterized by 
a single power of the logarithm ln(/i/(/)o) ~ l/<?. A part of them arises from the exchange of non-static magnetic and 
static electric gluons Another part is due to the contribution from the regular quark self-energy 0, 0] . Vertex 
corrections were reported not to contribute to subleading order |3 ■ The coefficient /3 in Eq. (JIJ determines the 
constant b in Eq. J2J). The third term in Eq. Q summarizes sub-subleading contributions of order 0(g 2 ) with neither 
a collinear nor a BCS logarithm. 

In principle, on the quasi-particle mass shell the gap parameter is an observable quantity, and thus must be 
independent of the choice of gauge. However, the mean-field approximation to the QCD gap equation may violate 
this requirement. It was arg ued in Refs. [EII1II3 that, in mean-field approximation, gauge-dependent terms enter 
the QCD gap equation at sub-subleading order. The authors of Ref. [21j confirmed that this is indeed the case for 
Coulomb gauge and when taking the gap parameter on the quasi-particle mass shell. In contrast to this result, it was 
pointed out [a. l22| that, in covariant gauge, the gauge-parameter dependence already shows up at subleading order. 
Denoting the gauge parameter in general covariant gauges by £, there is then an additional factor exp(3£/2) to the 
prcfactor b in Eq. @. 
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It is a priori not clear why it should depend on the choice of gauge, at which order gauge-parameter dependent 
terms appear in the QCD gap equation in mean-field approximation. A posteriori, this result is maybe not that 
puzzling at all, since the mean-field approximation corresponds to a resummation of a particular class of diagrams 
which may contribute at different orders for different choices of gauge. In any case, in order to remove the gauge- 
parameter dependence at subleading order in covariant gauge, one has to take into account corrections to the mean- 
field approximation. Gauge-parameter dependent terms will then only occur at sub-subleading order. The obvious 
correction so far not taken into account in the mean-field approximation is to replace bare qqg vertices by full vertices. 
It is the purpose of this paper to show how the inclusion of vertex corrections guarantees the gauge-parameter 
independence at subleading order in covariant gauge. 

As a matter of fact, it is not even necessary to compute the vertex corrections explicitly. We shall prove the 
gauge independence in a rather convenient way by making use of a Ward identity which relates the vertex to the 
inverse propagator. This approach has been frequently applied to show the gauge independence of physical collective 
excitations in thermal gauge theories, like hot QCD j2J, [24|. However, in a superconductor, the existence of a 
fermion-fermion condensate necessitates the use of the Nambu-Gor'kov (NG) basis to describe the propagation of 
quasi-particle excitations. It is therefore desirable to derive the Ward identity in the NG basis. In this paper, we 
derive this generalized Ward identity for color-superconducting quark matter and apply it to the gap equation, thereby 
showing that the gap parameter is gauge-independent to subleading order on the qua si-particle mass shell. 

We note that our approach is similar in spirit to that of Gerhold and Rebhan [2{|, who used generalized Nielsen 
identities to give a formal proof that the fermionic quasi-particle dispersion relations in a color superconductor are 
gauge independent, assuming that the 1PI part of the variation of the effective action with respect to the gauge 
parameter has no singularities coinciding with those of the quark propagator. Other Ward identities for color- 
superconducting quark matter have been derived by Miransky, Shovkovy, and Wijewardhana |2r| . but these identities 
were not suitable to see the gauge independence of the gap parameter to subleading order. 

In order to derive the generalized Ward identity, consider QCD with Nf quark flavors and quark chemical potential 
ji. The Lagrangian and the generating functional can be written in the NG basis, 



Z[J, H, H] = J [dA] [dtf] [d*] exp !j [£ + J°A% + + 



(4) 



The quark fields and their sources are given by 

*=(X)' * = $'^)' H=^y H = (rj,f} c ), (5) 

where the charge conjugate spinors ipc, "0c are defined through 4>c = Cip , ip G = -0 T C with the charge conjugation 
matrix C — ij 2 jo- The inverse free quark propagator (including the chemical potential fi) is 

where S$n = i^d^ + //70 — m, Sq% 2 = *7^/i — Wo — m - T ne c °l°r part of the qqg vertex is 

T a = ( Q _ T aT J ) (7) 

where T° are the generators of SU(3) C . The factor 1/2 in front of the first two terms in Eq. © accounts for the 
doubling of the quark degrees of freedom in the NG basis. The gluon field and the field strength tensor are denoted 
by A^ and F" , respectively. The last term in £ is the gauge fixing term, and J° is the source term for the gluon field. 

We denote the space-time integration as J x = J^ T dr J v d 3 x . We neglect ghost terms, as ghosts do not contribute 
to leading and subleading order in cold, dense quark matter Q. 

The generating functional Z is invariant under the following infinitesimal SU(3) C gauge transformation 

5* = id a T a ^ , 

5V = -i6 a ^T a , (8) 
5AI = ^d„e a + f abc Al6 c . 
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FIG. 1: The quark self-energy in the Dyson-Schwinger equation. The dashed line represents the hard-dense-loop (HDL) dressed 
propagator, the solid line is the quasi-quark propagator, and the black blob is the full qqg vertex. 



The variation of Z with respect to the infinitesimal parameter 9 a of SU(3) C gauge transformations vanishes, which 
yields the following identity: 



~ [d 2 dM)+gr bc [d,d„Dl»(X 7 X) + (Aftd^JQ] } - « +9f ahc J h M») +ig [HT a (V) (*)T»H] = , 

(9) 



where the average is taken in the functional sense 



(F(A, = ! J[dA][d*]l<M]F(A,*,9) ™p{J x [C + J^A* + H* + VH]} , (10) 

and = S 2 W/5J b 5J£ is the two-point Green's function, where W = In if is the generating functional for connected 
Green's functions. In the hard-dense-loop (HDL) approximation, the two-point function is symmetric in the 
color indices b and c. Consequently, the term proportional to Dj^ in Eq. vanishes, because the structure constants 
abc multiplying the two-point function are antisymmetric. In a color superconductor this does not need to be true 
but the violation of symmetry occurs for gluon energies and momenta of order ~ 6n. In the QCD gap equation, 
this range of gluon energy and momentum contributes only beyond subleading order pj . 

We can use the effective action T[A, = W - f x [J£A% +^H] to rewrite the identity J3J, i.e., we replace 

J£ by -ST/ 6 A*, ~H by ST/S^, and H by -5T/6l>. Note that, from now on, A", and * denote the expectation 
values of the gluon and quark fields. Taking the functional derivative S/S^B(Xi) from the right and S/S^(X2) from the 
left, and using the fact that the expectation values of the Grassmann- valued quark fields always vanish, \& = ^ = 0, 
we obtain the following identity 

_ 1 \QX S ab _ gf abc A c( X) l _ S 3 T[A,0,0] = tf(4) ( _ _6 2 T{A,0,0} 

5 Lm ' n S^{X 1 )SA b l (X)5^f(X 2 ) V 2 'S^{X 1 )S^(X 2 ) 

-S^X-X^J^S . (11) 
1 ' 8^{X l )5^{X 2 ) V 1 

The expectation value A ^ of the gluon field on the left-hand side is not necessarily zero in a color superconductor 
I27L l28j ] . However, in a two-flavor color superconductor, its value is of order 0q/ (g 2 /z) [2S|. i.e., much smaller than 
the derivative which is of the order of the gluon energy, po ~ cf>o, or even of the gluon momentum, p ~ (g 2 /i 2 ipo) 1 ^ s 
0. Therefore, in the following we may safely neglect the term proportional to A^ on the left-hand side of Eq. Ill IB . 
We now rewrite Eq. I jlljl in momentum space, 

P^(K,P,K + P)=g[T a S- 1 (K + P)-S- x (K)T a ] , (12) 

where T£(K, P,K + P) is the Fourier transform of 5 3 T /S^(Xi)8A® (X)5 i £(X 2 ), representing the full qqg vertex, while 
t S _1 (-ftr) is the Fourier transform of 5 2 T /S'^(Xi)S'i>(X 2 ), i.e., the full inverse quark propagator. We note that the 
above identity is a 2 x 2 matrix identity in the NG basis, and thus a generalized Ward identity. Equation l|12l) is 
rather similar to the Ward identity derived by Nambu to prove the gauge invariance of the Meissner effect in ordinary 
superconductors j^jji except that T a assumes the role of the Pauli matrix T3 in NG space. 
The inverse of the quark propagator S^ 1 in Eq. (|12f) is defined as 

S- 1 = So 1 + E , (13) 

where E is the quark self-energy in the NG basis. In momentum space, the Dyson-Schwinger equation for the quark 
self-energy reads 

m ) = -9 j n S(K + P) T b lu D%(P) . (14) 
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The right-hand side of this equation is shown diagrammatically in Fig. ^ Here, is the HDL-resummed gluon 
propagator and J p = T^2 n J d 3 p/(27r) 3 , where the sum runs over all fcrmionic Matsubara frequencies ui n = (2n + 
\)nT. The first vertex in the above equation is the full one and the second vertex is the bare one. 

We now show with the help of the generalized Ward identity (|12[) that, on the quasi-particle mass shell, the quark 
self-energy <|14ll does not depend on the gauge parameter £ in covariant gauge. To this end, consider the contribution 
from the gauge-dependent part of the gluon propagator to the quark self-energy, 

g£ j^S(K + P)T a 7v —^ ■ ( 15 ) 

Using the generalized Ward identity l|12|l . this becomes 

MK) = 5 2 C J [T a S-\K + P)- S-\K)T a ) S(K + P) T» 7 „^ . (16) 

Now we put the external momentum K on the quasi-particle mass shell, which is defined as 

5- 1 (^)* on _ shell = , * on _ sheU< S- 1 (Js:) = 0, (17) 

where V^on-shcii is the on-shell quasi-particle wave function. Then, when sandwiching Eq. I|16|) between on-shcll wave 
functions, the second term in brackets in Eq. IjlBjl vanishes. We arrive at 

MK)= 9 H\l / p ^5=0, (18) 

where 1 is the 2x2 unit matrix in the NG basis. The right-hand side vanishes because the integrand of the P-integral 
is an odd function of P. Hence, we see that the contribution from the gauge-dependent part of the self-energy Ijl4(l 
vanishes on the quasi-particle mass shell. Therefore, the gap parameter as well as the quasi-particle dispersion relation 
are gauge-independent in a color superconductor. 

This result holds strictly speaking only to subleading order in the QCD gap equation. The reason is that the 
generalized Ward identity l|12f> is not an exact identity, but only holds under the following approximations. The 
first approximation is that we neglected the contribution of ghosts. For QCD, which is a non-Abelian gauge theory, 
in general one has to take into account ghosts and use the BRST transformation and the Slavnov- Taylor identity 
or, equivalently, generalized Nielsen identities |U|2^|. However, in dense quark matter, diagrams involving ghosts 
can be neglected at least to subleading order in the QCD gap equation Q|. The second approximation is that we 
neglected the contribution from D^" in Eq. arguing that the gluon propagator is symmetric in the color indices. 
This is true for an HDL-dressed propagator, but not necessarily for the gluon propagator in a color superconductor 
Q. However, terms which violate the symmetry of should also be negligible up to subleading order in the gap 
equation. Finally, the third approximation was to neglect the expectation value of the gluon field in the generalized 
Ward identity l|12l) . In Ref. 28] it is shown that actually ^ in order ensure color neutrality. However, at least 
in a two-flavor color superconductor, A a ~ 0Q/(g 2 /i), i.e., the expectation value of the gluon field is negligible up to 
subleading order in the gap equation. 

In summary, we derived a generalized Ward identity from QCD for dense, color-superconducting quark matter. 
The identity implies that, on the quasi-particle mass shell, the gap function and the quasi-particle dispersion relation 
are independent of the gauge parameter in covariant gauge up to subleading order. There is one potential caveat 
to this statement. What we have shown is that, to subleading order, the gauge dependence of the quark self-energy 
(f 141) arising from the gauge- dependent part of the gluon propagator vanishes on the mass shell. In principle, however, 
other gauge-dependent terms arise from the gauge dependence of the full vertex T a ^ when combined with the physical 
part of the gluon propagator in Eq. I|14|) . In order to show that this is of sub-subleading order in the gap equation, 
one actually has to compute the gauge-dependent part of the vertex explicitly. This calculation can be simplified 
considerably with the help of the Ward identity l|12fl . for more details see Ref. [3(J ■ The result is that also this 
contribution is of sub-subleading order in the gap equation. 

Our result shows that in order to obtain a gauge-independent gap function up to subleading order, one has to use 
the full vertex as well as the full fermion propagator in the NG basis. A consequence is that the prefactor exp(3£/2) 
to the gap function found in the mean-field approximation 0, [2^ will be removed by contributions from the full 
qqg vertex when taking the gap function on the quasi-particle mass shell. An explicit diagrammatic proof of this 
statement will be presented elsewhere |3Cj . 
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